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Abstract 

The use of the Boson Loop Expansion is proposed for investigating the static prop- 
erties of nuclear matter. We explicitly consider a schematic dynamical model in 
which nucleons interact with the scalar-isoscalar a meson. The suggested approx- 
imation scheme is examined in detail at the mean field level and at the one- and 
two-loop orders. The relevant formulas are provided to derive the binding energy 
per nucleon, the pressure and the compressibility of nuclear matter. Numerical re- 
sults of the binding energy at the one-loop order are presented for Walecka's a-u> 
model in order to discuss the degree of convergence of the Boson Loop Expansion. 
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1 Introduction 



Functional integral methods have been extensively used over the years to 
provide useful approximation schemes in many different areas of physics, in 
both the perturbative and the non-perturbative regimes [1-6]. The nuclear 
many-body problem has been widely investigated starting from renormaliz- 
able, relativistic Lagrangians based on the so-called quantum hadrodynamics 
(QHD) [7]: calculation schemes have been explored, which go beyond the tree- 
level approximation and thus involve the dynamics of the quantum vacuum. 
Being a strong-coupling model, no obvious asymptotic expansion exists for 
QHD. Hence it is not guaranteed that a reliable, convergent expansion, al- 
lowing for systematic refinements of the theoretical predictions, exists for this 
model. 
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Some approaches are based on the so-called Loop Expansion, which can be 
derived from the exact path integral formulation of QHD [4]: it is obtained 
by expanding the QHD effective action around its classical value and by con- 
sidering terms with increasing number of quantum loops [1,2]. Formally, the 
resulting Loop Expansion is thus a power series in h, though Planck's constant 
is merely a bookkeeping parameter [2]. One can also consider it as a perturba- 
tive expansion in powers of the coupling constants g a , g^, etc, but keeping in 
mind that it is non-perturbative in the mean fields, since it involves Hartree- 
dressed nucleon propagators. The one-loop approximation of the Loop Expan- 
sion corresponds to the Relativistic Hartree Approximation [8]. In Refs. [9] full 
two-loop corrections were evaluated for QHD-I (i.e., Walecka's a-cu model) and 
criteria for "strong" and "weak" convergence were discussed. Two-loop terms 
involving the vector meson oj were found to be too attractive for standard 
values of the coupling constant g u , thus making two-loop corrections to the 
nuclear matter binding energy too large with respect to the well known mean 
field results. This result precluded strong convergence. After tuning the QHD- 
I parameters, the correct saturation point of nuclear matter at the two-loop 
order was obtained at the price of drastically reducing g w . At the same time 
g a was significantly increased, the dynamics of the system being led almost 
entirely by the scalar a meson. It was thus concluded that, at two-loop order, 
the Loop Expansion is not even weakly convergent. 

To suppress the unwanted high-momentum structure of the hadronic vac- 
uum loops, short-range correlations and vertex corrections have also been 
considered in calculations at the two-loop order [10-12]. These investigations 
showed a clear improvement over the results of Ref. [9], which were carried 
out with structureless nucleons. Although obtained with a "reasonable" fit of 
the coupling constants, the results turned out to have a strong dependence on 
the model, not phenomenologically constrained, adopted for the short-range 
correlations or the vertex corrections. Moreover, the renormalizability of the 
theory was lost. 

The above failure of the Loop Expansion could be due to its inability in de- 
scribing the quantum vacuum dynamics. Indeed it is derived as a perturbative 
expansion in the "large" QHD-I couplings, while "small" couplings are re- 
quired to assure some convergence of the Loop Expansion predictions [13]. 
One could also claim that the QHD model itself, based on point-like hadrons, 
is not appropriate to include quantum vacuum effects. 

In order to test the validity of these hypotheses, alternative loop expansions 
based on QHD, non-perturbative in both the couplings and the mean fields, 
should be studied. In this context we mention here the so-called Modified 
Loop Expansion or Boson Loop Expansion, which is an example of a non- 
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perturbative approach and in addition preserves renormalizability 1 . 

The Modified Loop Expansion, first introduced by Weiss [15], has been later 
employed in Ref. [16] together with a chirally symmetric linear a model for 
nuclear matter calculations. In this approach, fermionic degrees of freedom are 
integrated out in the generating functional. This leads to an effective action 
in which meson loops are included order by order in H and, at each order h n in 
the meson loops, baryon loops are summed to all orders in h. The zeroth-order 
of this expansion corresponds to the Relativistic Hartree Approximation, the 
modified one-loop level to the Relativistic Random Phase Approximation. 
Hence the first-order Modified Loop Expansion corresponds to the two-loop 
level of the Loop Expansion, but with meson propagators dressed in Random 
Phase Approximation. We are not aware of detailed calculations applying this 
scheme to QHD-I. 

The so-called Boson Loop Expansion (BLE), which we shall consider in the 
present paper, has been proposed as a suitable tool to formulate consistent 
approximation schemes for non-perturbative calculations in different issues of 
nuclear physics [17-23]. The underlying idea, as in MLE, amounts to integrate 
out the nucleonic degrees of freedom within a path-integral formulation for the 
dynamics of interacting nucleons and mesons. The loop expansion is then car- 
ried out on the resulting bosonic effective action. This approach has been ap- 
plied to the calculation of the electromagnetic nuclear response functions and 
provided a suitable framework to understand the reaction mechanisms [21]. 
It has also been successfully used for the evaluation of the hypernuclear weak 
decay rates [23] . 

Actually, while the BLE has the merit of rigorously preserving sum rules and 
general field theory theorems, the underlying dynamics still needs to be more 
firmly settled. We intend to apply this method to the investigation of other 
observables, such as the static properties of nuclear matter and neutron stars, 
topics of considerable interest and yet under debate after decades of valuable 
works. It is worth noticing that this approach is particularly suited to deal with 
a relativistic description and, in addition, it is stable against sizable increases 
of the nuclear density. 

Particular attention has to be payed to the adopted dynamical model. In a 
relativistic frame, QHD models are often employed at face to calculations 
based on realistic nucleon-nucleon potentials like the Bonn meson exchange 
model [24]. These two schemes started from quite different points of view 
and originally provided sizably different results: however in recent times the 

1 We recall that another approach has been recently suggested to overcome the dif- 
ficulties inherent with the description of the vacuum dynamics in QHD: it amounts 
to utilize QHD as a non-renormalizable, Effective Field Theory in calculations at 
the tree-level [14]. 
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evolution of calculations in both approaches has significantly reduced the gap. 
QHD is mainly based on the a-ou interference, while the Bonn potential, in 
addition to the exchange of several physical mesons, replaces the a meson 
exchange with a two-pion exchange together with the simultaneous excitation 
of one or two intermediate nucleons to A resonances (box diagrams). From a 
phenomenological point of view, this seems to favor the Bonn approach for the 
calculation of the equation of state of nuclear matter. However, a G-matrix 
calculation based on the Bonn potential is quite cumbersome and some self- 
consistency effect of relativistic origin is generally neglected. In this context, 
calculations are usually limited to the Lowest Order Bruckner Theory, which 
implies lowest order in the nuclear density; this limitation is acceptable in 
ordinary nuclear matter calculations, but less so in addressing dense baryonic 
matter such as in neutron stars. 

In order to clarify the fundamental problems underlying the BLE (namely 
the issue concerning the renormalization procedure), in this paper we develop 
the formalism within a simple dynamical model containing nucleons and a 
mesons only. We shall provide the relevant formulas which are needed to derive 
the static properties of nuclear matter up to the two-boson-loop level. In 
view of a future, comprehensive application to QHD, as well as to give a 
first insight on the convergence properties of the proposed expansion, here 
we present numerical results for the nuclear matter binding energy at the 
one-loop order when adopting a g-uj model. The complete BLE formalism for 
QHD-I, together with systematic numerical calculations up to the two-loop 
order, will be presented in a forthcoming paper, also devoted to discuss the 
phenomenological aspects of the present approach. 

The paper is organized as follows. Section 2 contains a brief introduction to the 
general formalism employed in the functional approach for the evaluation of 
the static properties of nuclear matter. In Section 3 a bosonic effective action is 
derived for a system of interacting nucleons and a mesons. The "elementary" 
vertices entering into this effective action are analyzed in some detail in Section 
4. In Section 5 the renormalization procedure for the bosonic action is carried 
out. In Section 6 we deduce the BLE as a Semiclassical expansion around 
the (mean field) solution given by the stationary phase approximation. Some 
numerical results for the a-oo model are presented in Section 7. Finally, in 
Section 8 we comment on the future perspectives opened by the present work. 



2 The functional method - General formalism 

To construct the formalism we consider a simplified situation where a scalar- 
isoscalar meson (the a) interacts with a massive fermion field (the nucleon). 
The introduction of the u meson, to complete the standard scheme of QHD 
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induces a mixing with the a meson and a problem of renormalizability, which 
can be solved by means of the so-called Stueckelberg gauge [25] and not simply 
invoking the current conservation, as usually stated [7]; this will be discussed 
in a forthcoming paper. We shall see in the following that already the "simple" 
scalar field needs a careful renormalization procedure. 

The most general Lagrangian density reads 

— 1 1 — 

C = - m)ip + -d^a d»o - -mla 2 + gipaip - U(a) 

= J d%^(x)S» 1 (x-y)^(y) + ±jd 4 ya(x)D 1 (x-y)a(y) W 
+ gip(x)a(x)ip(x) - U[a{x)\ , 

the two expressions differing by an irrelevant four-divergence. In the above 

S \x-y) = ( i @ x -m)6\x-y) , (2) 
D \x -y) = -{U x + ml) 5\x - y) , (3) 

formally define the free nucleon and a propagators, So and D respectively, in 
coordinate space, and 



t^) = ^ 3 + ^ 4 = ^- 3 + ^V. (4) 



In momentum space the nucleon propagator in nuclear matter with Fermi 
momentum kp and the free a propagator read, respectively: 

It + m j 9(\k\-k F ) 6(k F -\k\) 1 \ 

° l ' 2E k [ko-Ek + ir] k -E k -ir) k + E k - t V J ' {) 

DM = \— - , (6) 

q z — m l a + le 



where E k = Vk 2 + m 2 . It is also useful to introduce a symbol for the free 
fermion propagator in the vacuum, S 0v = S \ kF=0 , thus splitting Eq. (5) as 
follows 



S (k) = S 0v (k) + S 0M (k) , (7) 

Sov(k) = - (8) 

f — m + ir\ 

S 0M (k) = (It + m)—6(k F - \k\)5(k - E k ) , (9) 
the latter being the density-dependent part. 
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The classical action corresponding to the Lagrangian (1) is given by: 



where the a and nucleon fields have to be interpreted as scalar and Grass- 
mann variables respectively. In a path integral approach the whole dynamics 
is deduced from the generating functional 



where Af is a normalization constant while 77, fj and J denote the external 
(classical) sources of the nucleon and a meson fields. From Z one derives 
the nucleon Green's functions by functional differentiations with respect to 
i] and fj. The a Green's function is obtained by derivatives with respect to 
J. Other external sources coupled to composite fields could also be added: 
for instance, couplings of the kind ?p(x)il)(x)p(x) permit the evaluation of the 
response function through a double derivative of Z with respect to p. 

The static properties of the system follow from the partition function Z((3) 
which can be obtained from the generating functional of Eq. (11) by: 

(1) setting the external sources to 0; 

(2) performing a Wick rotation, namely replacing, in the exponent, the time 
integration (extended from —00 to 00) with an integration over the imag- 
inary time r = it in the interval [0,/3] with (3 = H/ftsT) = 1/T; 

(3) replacing the Hamiltonian H = J d 3 x7i(x) with H — pN, p being the nu- 
cleon chemical potential and N = J d 3 x , ijj(x)'j ilj(x) the nucleon number 
operator; 

(4) setting M = 1. 

By applying these rules, one gets 



The Euclidean Lagrangian Ce is obtained by replacing t with — ir in the 
Minkowskian Lagrangian (1) and reads 



where we have denoted the Euclidean space-time coordinates with an upper 




(10) 





(12) 
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bar [x = (r, x), y = (r y ,y)] and TC is the (Euclidean) Hamiltonian density. In 
the zero temperature limit Sq^ and Dq^ are given by 



S^{x -V) = -(l°d Tx + 1 7 • V, + m) 5\x - y) , (14) 
Doe(x -V) = -(-d 2 Tx + ml) 5\x - y) . (15) 

The Euclidean propagators in momentum space follow from Eqs. (5), (6) with 
the replacements k — > ik + ji and q — > iq and take the form 



S 0E (k) = , (16) 

7o(ik o + a0 - 7-k-m 

D 0E (q) = - 2— 4t 2 • ( 17 ) 



Once the partition function is known, the ground state (T = 0) energy 
of nuclear matter follows from the relation 

E - W = ("! ^ + = >to (-^ + M) . (18) 

A = (TV) being the total number of nucleons. Note that the nucleon chemical 
potential is related to the Fermi momentum by 



lim n = E* F = Jk% + (m*) 2 , (19) 

/3— >oo v 

where m* is the nucleon effective mass at the mean field level [see Eq. (38)] and 
follows from the dynamics underlying the Lagrangian (1). In the next Section 
we shall introduce an approximate scheme to evaluate the ground state energy 
and therefore the binding energy per nucleon: 

B.E. E e 

— : — = —r — m = m , (20) 

A A p 

where e = E/Q and p = A/Q = 2kp/(3n 2 ) are the energy density and the 
number density of nuclear matter, respectively. Note that we are implicitly 
considering symmetric nuclear matter (/cp n = kp p = kp) but all our results 
can be easily adapted to the asymmetric case. 

Once e is known, other observables are accessible. Pressure (and hence the 
equation of state) and compressibility are indeed given by the thermodynam- 
ical relations: 

„ 9 d ( e\ de , . 
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and 




(22) 



Other quantities, such as for instance the symmetry energy, could also be ad- 
dressed, but they involve the dynamics of isovector mesons, which has been 
recently studied in Refs. [26,27], and will be considered by us in future inves- 
tigations. 

Within the present path integral method, the problem of evaluating the static 
properties of nuclear matter requires the elaboration of approximation tech- 
niques in order to compute the functional integral which defines the partition 
function (12). 



3 The Bosonic Effective Action 

As stated in the Introduction, we shall concentrate here on the approximation 
scheme given by the Boson Loop Expansion [17, 21]. As a first step to derive it, 
one has to introduce a bosonic effective action corresponding to the Lagrangian 
density (1) but no longer displaying fermionic degrees of freedom. This is the 
specific purpose of the present Section. 

At variance with Ref . [8] , we present the theory in a heuristic way, disregarding 
in the present Section the renormalization problem and the a self-interaction 
terms — which in turn add non-trivial complications to the renormalization 
problem — and leaving for a second step (Sec. 6) the details of the renormal- 
ization procedure; the latter, indeed, affects the binding energy in the medium 
already at the mean field level [8]. We also notice that the classification of the 
counterterms in the usual Loop Expansion and in the BLE is different. 

Since we shall focus our attention only on the boson-like observables of Eqs. (18)- 
(22), we can safely put f) = i] = in the generating functional of Eq. (11). 
The functional integral is then converted into an integral over bosonic variables 
only, by explicitly integrating out the fermionic field; one obtains 



By using well known properties of the following Gaussian integral over Grass- 
mann variables: 




(23) 




the bosonic effective action A^ s entering into Eq. (23) turns out to be 
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= \J d 4 xd 4 ya(x)D 1 (x - y)a(y) - zTr \og[-i(S^ + ga)\ (24) 
= - J d 4 x d 4 y a{x)D \x - y)a(y) + A* + *Tr £ L-i- [S o9 a] n 



n=l 
9 



\ J d 4 xd 4 ya(x)D^(x - y)a{y) + A* - iTr J 



d\ Act 



X So 1 + Xa ' 
u 



The term 



A* = -iTrlogl-iSo 1 ] (25) 
is related to the functional generator of non-interacting fermions by 

e iA$ = j JJ^p^^ij^xdS^Sg^x-y^iy) _ ^6) 

Let us stress that, although the fermionic degrees of freedom have been inte- 
grated out, their influence on the generating functional (23) remains unaltered. 

Note that in (24) space-time integrals are also contained in the traces, in 
addition to a sum over spin and isospin. The trace of the logarithm acquires a 
meaning only through its series expansion (second line of Eq. (24)). Moreover 
Tr[S ga] n is a short notation for 

Tr[S ga] n = Ag n J d 4 x x . . . d 4 x n ^\x u . . . , x n )a(x 1 )a(x 2 ) . . . a(x n ) , (27) 

where we have defined 

^\xi, ...,z n ) = S (xi - x 2 )S (x 2 -x 3 ) ... S (x n - xi) , (28) 
and the factor 4 originates from the spin-isospin trace. 

It will be also useful to introduce the vertices evaluated in the vac- 
uum, denoted by ^3 ™^: they are obtained from Eq. (28) with the replace- 
ment So — > S 0v , while the corresponding density-dependent contributions are 
m( n ) _ m( n ) <n( n ) 

+>0M — +>0 +>(k> • 

Their Fourier transforms are also needed. Here the independent variables are 
n — 1 and we get 



?0 n) (Pl,P2, 

d 4 p 



■Pn-l) 



(29) 



/dp 
j^^S (p)S (p + pi)S (p + pi+p 2 )... S (p + pi H h Pn-l) ■ 
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The Euclidean version of the bosonic effective action is 



A* eS ([3) = lHog[-S -i] + A B Ee M , (30) 
A B Ee M = 5 / ^rr/Vi.r)M, / M.r - j/)cr(j/) - Tr £ fcH" 

27 n=l n 

= \ J d*xd*ya(x)D^(x - y)a(y) + Tr J f \a ' 

and the partition function takes the form 

Z((3) = Z ((3) j D[a}e A *^ . (31) 

Here, 

Z (p) = deff-S^ 1 ] = f D$, il>\ e Jd 4 xd^mSo^-y)m ( 32 ) 

is the partition function for an assembly of non-interacting nucleons and the 
remaining functional integral in Eq. (31) represents a system of self-interacting 
bosons. 

The interaction terms of Eq. (27) are built up with closed fermionic loops, 
either in the Minkowskian or Euclidean space, and play a central role in the 
present treatment. It has been shown in Refs. [18,28] that, within a non- 
relativistic kinematics, they can be either evaluated explicitly or at least re- 
duced to the evaluation of a one-dimensional integral. Note however that the 
difficulties met in Refs. [18, 28] mainly derived from the kinematical singu- 
larities in Eq. (27). The use of the Euclidean metric shifts all of them in 
the complex plane, considerably simplifying the numerical integrations. Since 
these interaction terms represent the building blocks of the BLE, they require 
a careful treatment. Relativistic calculations are much more sensitive to the 
dynamics of the mesons with respect to non-relativistic case. In the latter, 
indeed, at least in most practical cases, the mesonic dynamics is frozen (the 
meson exchange is reduced to a static potential) and it is reflected in some 
spin operators which are easily traced out. On the contrary, in a relativistic 
approach convective currents cannot be neglected without violating Lorentz 
covariance. We shall discuss in the next Section how some archetypal cases 
may be constructed and how one can handle more complicated dynamics. 

To evaluate the generating functional (23) or the partition function (31), we 
adopt the semiclassical expansion. Its lowest order, i.e., the mean field level, 
corresponds to the stationary phase approximation (saddle point approxima- 
tion in the Euclidean space). At this stage, one requires the bosonic action to 
be stationary with respect to small variations of the field a: 

s jm =0 . (33 ) 

0(7 
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The following equation of motion is thus obtained: 



J d 4 yD 1 (x-y)a(y) = -*Tr £ (-l) n [5 H n_1 ^ = ^ Tr 1 + g Qg(7 • ( 34 ) 

One recognizes that cr cannot be vanishing at the mean field level, since the 
first term of the sum in Eq. (34) is just a fermionic line closed onto itself, 
namely the tadpole 

igTrS = 4ig J d 4 y S (y - y) . 

The above quantity is obviously divergent but the usual renormalization tech- 
niques, described in Sec. 5, make it finite and non-trivial [8,29]. This occur- 
rence emphasizes the peculiar role played by the scalar-isoscalar field, which 
is indeed the only one generating a tadpole, even in the vacuum, where any 
other meson field has vanishing expectation value. 

For static, infinite nuclear matter, the mean field solution is uniform in space 
and time, thus the lhs of the field equation (34) reduces to 

J d 4 yD^(x - y)a{y) = -(□* + ml) a(x) = -m 2 a a . (35) 

In the "no-sea" approximation, the mean field a is then obtained by solving 
the following self-consistency equation (we introduce the scalar density p s = 
{ipip} and its density-dependent part p s u = {^i])) m)'- 

v = ~^PsM = -i—z-Tr S H m = -4i— - / — ,, S H M(k) 
ml ml ml J (2vr) 4 

9 m* ( k F + E* F \ 1 ] 

= — — k F E F - (m log — , 

m z y m* J 

where Shm is the density-dependent part of the Hartree-dressed nucleon prop- 
agator 

Sh = , (37) 

b + go- 
that is obtained from the definition of Sq in Eq. (5) after replacing m with the 
nucleon effective mass: 

m* = m — ga . (38) 



From a diagrammatic point of view, self-consistency amounts to account for 
diagrams like those of Figure 1, in which each nucleon propagator is dressed 
by the Hartree self-energy. 

In view of an application of the present scheme to QHD, one has to consider 
the contribution of the vector-isoscalar meson uj as well. The ui mean field 
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Fig. 1. Mean field level nucleoli self-energy diagrams for a exchange. 

is obtained as uj^ = (^ w / TO w)/°^o ? without any self-consistency condition 
because the nuclear density 



P = (^70-0) 



3tt 2 



is now implied instead of the scalar density p s . The relativistic (Lorentz con- 
traction) effect stemming from the difference between the two mean values 
was already noted long time ago by Lee and Wick [13] and ultimately entails 
a restoration of the chiral symmetry at sufficiently high density in the Wigner 
mode. 



Since Eq. (36) provides a non-vanishing solution for the mean field a, in 
order to apply the standard technique of the Semiclassical approximation one 
needs to define a new scalar-isoscalar field a\x) = a(x) — a with vanishing 
expectation value at the mean field level, a' = 0. By rewriting the generating 
functional (23) with the above change of variable and then renaming a' as a, 
the bosonic effective action (24) becomes 

AM = ~^D-'a - -f{2ao + a 2 ) + A* + *Tr £ { -^-[S g(a + a)]» (39) 

22 n=l H 

and the coupling to the external field changes to (a + a) J. 2 

2 Note that a compact notation has been used for the first two terms of the effective 
action (39). Their explicit expression reads: 

\ J d 4 xd 4 y [a(x)D 1 (x-y)a(y)-m 2 j\x-y)(2a(x)a + a 2 )] . 
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Next, one can write 



tTr ^ K -^[S g(a + a)] n 

n=l H 

00/ I \n 00 / 1 \n 

= ,Tr5:^^[^a] n + ^r5:^^[^ar, (40) 

as it can be easily verified keeping in mind Eq. (37) and by rewriting the series 
in terms of the corresponding log's. By using Eq. (36), the n = 1 term in the 
first sum of the rhs of Eq. (40) reads —igTrSnv = m 2 a aa , so that no linear 
term in a enters the effective action: 

00 (_l)n 

+ A B + 1T1 ]T K —^[S H ga] n - iTi log[l + S ga] . (41) 
„= 2 n 

Let us remark that the constant terms in are irrelevant in evaluating the 
generating functional, but they are significant in the partition function. By 
using the definition (25), one can recognize the following reordering 

A B - *Tr log[l + S ga] = -tTr logf-^ 1 ] . (42) 



We then consider the partition function. The part of the Euclidean action orig- 
inating from the constant term (42) now replaces the factor Zq((3) of Eq. (32). 
The difference between Eq. (42) and Eq. (25) lies in the replacement m — > m*, 
namely Soe —> She = (Sq E + go-)' 1 , She being the Euclidean version of Sh- 
Therefore we can write 

Z(P) = Z* (P)jDW]e^^ m , (43) 

with 

Z*(p) = det[-S HE ] = J DfcM] e^^^^H^-y)^) , (44) 

and 

A B Ee M = ^D^a - Qp^a 2 - Tr £ [ —^- [S HE ga\ n . (45) 

2 2 n=2 71 

In this expression, space-time integrations are understood for the free a con- 
tribution, but in the second, constant term in the rhs, the time {(3) and volume 
(Q) integrations have been made explicit. 
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4 The elementary vertices of the bosonic action 



We analyze now the elementary vertices corresponding to Ti[S H go] n and enter- 
ing the bosonic effective action (41). When embedded in a Feynman diagram, 
a term Tr [S^gcx]™ is represented by a fermionic loop with n external points; a 
boson propagator can be attached to each one of these points, including the 
spin and isospin matrices which pertain to the specific boson. Isospin traces 
only factorize out a real number, but the possible presence of 7 matrices, as 
it is the case for the u exchange, entails a non-trivial dependence upon the 
momenta, which is usually neglected in non-relativistic calculations. In this 
section we only consider some archetypal structures. More realistic situations 
need to be dealt with individually. 

Let us first fix the kinematics. A loop with N + 1 incoming or outgoing mesons 
will depend upon N momenta pk, while pn+i = — J2k=iPk is fixed by mo- 
mentum conservation. To simplify the notation, it is convenient to define N 
auxiliary momenta according to 



Qo = , 



Vi = J2Pi (i = l,-..,N) . 

3=1 



(46) 




Fig. 2. The kinematics of the generic interaction term in the expan- 
sion of the bosonic effective action (41). 

The kinematics is clarified in Figure 2. With these definitions the generic 
vertex in momentum space reads: 

/d^p 
j2^4°o(P' Vo)Sh(p + qo)Oi(p, qi)S H (p + q x ) . . . N (p, q N )S H (p + ftv) , 

where the Oi are matrices whose particular form depends on the bosonic 
field(s) one considers. 

Actually, even the simplest case, namely the exchange of scalar-isoscalar 
mesons, is rather involved, because the traces over the 7 matrices are non- 
trivial. 
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In order to handle this case, we introduce a set of auxiliary functions as fol- 
lows. First we rewrite the nucleon propagator in nuclear matter with Fermi 
momentum k F as 

S (k) = (ft + m)e (k) , 

where 

2E k \k - E k + ir] k - E k - ir] k + E k - ir] 
Then we introduce the auxiliary functions 

^ N+1 \q , ■ ■ ■ , q N ) = -i [ ^ ft ©o(p + Qi) , (48) 

which, at variance with the of Eq. (28), do not contain the matrix struc- 
ture of the propagators So- 

Following Refs. [18, 28] we integrate over p by closing the integration path 
with a half circle in the lower half plane, getting 




j=Q J [ZTT) i=0 Ah p+q .h p+qi 

0(\p + <h\-k F ) 
Qio — Qjo — \E p + qi — E p+q .) 
9(k F - Ip + q^l) 1 1 



X 



Qio - Qjo - (E p+Qi - E p+q .) - it] q i0 - q j0 + E p+q . + E p+q . - ir] J 

The rhs displays a sum of products. Each term in the product is in turn a sum 
of three pieces, with poles at q i0 — qj — (E p+qi — E p+q .) = (first term), q i0 — 
q j0 -(E p+q .-E p+q .)-ir) = (second term) and q i0 -q j0 +E p+q .+E p+q .-ir} = 
(third term). In the first term, at variance with the others, the singularity is 
removable and hence no prescription on how to handle it is needed. Indeed 
this singularity arises if we consider the i-th term of the product in the j-th 
term of the sum; but, by exchanging % and j, a term with the same singularity 
and opposite sign appears, thus removing the singularity. Being irrelevant, we 
can formally add to the denominator of the first term an infinitesimal factor 
sign(q iQ — qj Q ). In the second term the 9 functions impose E p+Q .—E p+q . < 
and the denominator can vanish only when qio~qjo < 0. This allows to replace 
the —ir] by sign(q i0 — qj Q ) and to combine together the first two terms, 
the 8 functions summing up to 1. Finally, the same factor can be ascribed 
a fortiori to the imaginary part of the last term, thus yielding the compact 
expression: 
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N 
3=0 • 



n'~«»( 9 „,..., 9N ) = -E/(g5n- (|p + <1,l " M 



2£, 



(50) 



1 



X- 



Qrn ~ Qjo + E p+qj + irj sign(g i0 - q j0 ) 



E 2 



If we are only interested to the in-medium contribution of the above ver- 
tices we have to subtract from (50) its value in the vacuum; this leads to the 
following expression, where no residual divergence is left: 



N 



n^)( go ,...,^) = E 

3=0 



(f'p * 0{k,, 



Ip + QjI) 



(27T; i=0 



2E. 



(51) 



p+qj 
1 



X- 



Qio - Qjo + E p+q 3 + iv sign(g i0 - <ljo) 



E P+qr 



The transition to the Euclidean space is even simpler, because in the replace- 
ment q i0 — > iq i0 +/i the chemical potential cancels out, while the +irj sign(g i0 — 
qjo) becomes irrelevant. One gets: 



n 



(N+l) 



N 



(So, 



/ d 3 P A @(kp - Ip + q^l) 



2E. 



(52) 



X- 



iQio ~ iQjo + Ep+qj 



/?2 



Relativity forces in the elementary vertices of the bosonic effective action 
many complications with respect to the non-relativistic case discussed in 
Refs. [18,28], and no general formula can be given in the relativistic case. 
In Ref. [30] a detailed study of the two-point vertices has been carried out. 
Here we shall illustrate with a simple, but yet non-trivial, example how the 
three-point vertex can be reduced to the evaluation of auxiliary functions 
n( Ar+1 )(g , • • • , Sat). Let us consider the emission of three a mesons from a 
three-point fermionic loop. The corresponding vertex is 



^ (3 W) = -*/(0i (53) 

Tr(^ + to) {jf> + j^i + to) (j> + ^ 2 + m) 
[p 2 — to 2 ± irj][(p + gi) 2 — to 2 ± ir)][(p + q?) 2 — to 2 ± irj] 

(the sign of the irj follows from the above discussion). By performing the traces, 
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the numerator of the integrand becomes 



4m 3 + Ylmp 2 + 8mp ■ qi + 8mp ■ q 2 + Amqi ■ q 2 , 



and with some algebra can be recast in the form 



4m{(4m 2 - q\ — q\ + q x ■ q 2 ) 



+ (p 2 - m 2 ) + [(p + qi ) 2 - m 2 } + [{p + q 2 ) 2 - m 2 ]} . 



By inserting this result in Eq. (53) we then get 



^ (3) ( gi ,g 2 )=4 m (4m 2 - q\ - q\ + q x ■ g 2 )n (3) (q u q 2 ) 

+4mn (2) (g 2 - qi) + 4mn (2) (g!) + 4mn (2) (g 2 ) . 



(54) 



Hence in this particular case a realistic three-point vertex is reduced to the 
evaluation of the auxiliary vertices II^ 3 '' and IT^. The procedure can be sys- 
tematically generalized. 



5 Renormalization of the Bosonic Effective Action 

The derivation of Section 3 only has a formal meaning since no renormalization 
procedure has been applied up to now. To make the above derived quantities 
meaningful we need 

(1) to add to the Lagrangian a set of counterterms and 

(2) to fix, by altering the Feynman rules, a regularization procedure. 

In the present case it is convenient to rewrite the Lagrangian (1) for bare 
quantities, denoted by the index 0: 



Bare fields and constants are linked to the renormalized ones by the relations: 



C = ^ (ift-m )^ + -d^a d^a --ml aQ+goiJ (j ij 





_ 3 

^30 — 2 A3 



ip = Z 2 2 ip 



1 



(56a) 

(56b) 
(56c) 
(56d) 

(56e) 

(56f) 



A40 — Z 5 Z 3 A4 . 
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By substituting we rewrite the Lagrangian in the form 




Co + d + C, 



(57) 



where 




-\5mW + (Z S - 1) {^rffVr - ^m> 2 } 
+(Z! - 1)#^ - (Z 4 - l)^x 3 - (Z 5 - 1)^<7 4 . 



(58b) 
(58c) 



(58a) 



According to the general theory of renormalization [4,31], each order of the 
perturbative expansion remains now finite. It is convenient for our purposes 
to fix only one expansion parameter, the coupling g; this can be achieved by 
putting A3 = m^a^g 3 and A4 = a^g* and attributing to 5m, 5m 2 a , (Z 2 — 1) and 
(Z 3 — 1) the order g 2 , to [Z\ — 1) and (Z4 — 1) the order g 3 and to (Z 5 — 1) 
the order g 4 . 

Let us now fix a regularization procedure: its choice in many-body calculations 
is largely immaterial, but for the present purposes and for practical calcula- 
tions it will be useful to fix a scheme. We shall simply adopt a cutoff scheme, 
according to 



(the square is just the minimal required power to force convergence in any of 
the diagrams occurring in the theory). 

Pauli-Villars regularization [4] will also be kept in mind, since it will help us in 
the power counting inside the fermion loops. The dimensional regularization 
adopted in Ref. [8], instead, appears less natural in the present many-body 
context; indeed the nuclear medium breaks Lorentz covariance and it becomes 
unclear to which dimension (3-dim space or 1-dim time) the regularization 
has to be applied. 

A short discussion of a few possible renormalization schemes, although well 
known [4], is required in this context. 

(1) One can choose the parameters in Eqs. (58a)-(58g) as the physical masses 
and coupling constants, so that at each perturbative order a counterterm 




(59) 
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(infinite or finite) must be introduced to keep them fixed (hard renormal- 
ization). This entails that the renormalization point must be on-mass- 
shell for any external line in evaluating the relevant diagrams. 

(2) One can instead remove from any elementarily divergent diagram its di- 
vergent part (soft renormalization). Dimensional regularization is usually 
adopted within this scheme. 

(3) Intermediate choices are also possible. For instance Ref. [8] chooses the 
renormalization point at q 2 = (thus generating some inconsistency in 
the definition of the a mass). 

(4) Finally, as suggested in Ref. [29], one can adopt as counterterms the val- 
ues of the elementarily divergent diagrams in the vacuum, considered as 
functions of the incoming and outgoing momenta. A rigorous treatment 
of the renormalization is then lost, since the ill-definition of the Green's 
functions is no longer restricted to a point, but the resulting scheme is 
finite. Moreover one disentangles, at least partially, the physics of the nu- 
cleon, lying in the form of the vertex functions in the vacuum, from the 
nuclear physics. When this scheme, referred to as "vacuum subtraction" 
scheme in the following, is applied to QHD, one of the main advantages 
is to avoid meaningless contributions from the vacuum. One should in- 
deed keep in mind that QHD is intrinsically an effective theory, with 
parameters ruled by nuclear dynamics. 

A detailed discussion of the elementarily divergent diagrams is now in order. In 
fact the procedure does not follow the usual path of perturbation theory, since 
loops containing or not a boson line play different roles and the requirement 
of being one-particle irreducible (1PI) looses its validity in the context of the 
BLE (in fact it refers now only to boson lines). Furthermore the tadpoles 
(referred to in the following as Hartree contributions) play a special role. 

We consider first the Hartree contribution to the nucleon self-energy of Fig. 3. 
It plays a peculiar role in the BLE, since it will eventually disappear after the 



Fig. 3. a) The Hartree contribution to the nucleon mass renormalization. b) and c) 
The elementarily divergent diagrams with one internal boson line. 

self-consistency requirement is forced; thus it does not increase the bosonic 
loop number, neither the fermionic one. It reads 





a) 



b) 



c) 



^Hv — o 




(60) 
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where p sv = {ijjip)\k F =o is the scalar density in the vacuum. We have em- 
phasized its dependence upon m since, as already heuristically anticipated in 
Section 3, the replacement of m with the effective mass is expected. The ex- 
plicit expression of p sv (m) within the regularization scheme given by Eq. (59) 
is provided in the Appendix. 

Next we consider the divergent diagrams with one bosonic loop. They are 
displayed in Fig. 3 and represent the Fock term of the self-energy (the name 
being borrowed from the standard many-body language) and the aNN vertex 
correction. The diagrams a) and b) of Fig. 3 remove the divergences from the 
nucleon self-energy. 

In full generality, in the vacuum the inverse propagator of the nucleon can be 
written in the form (we use the standard notation for the 1PI vertices and the 
corresponding generator, T) 3 : 

5 2 T 

r^(p) = jrir- = Z 2 {/>-m) -8m-T, v {p) , (61) 

where the self-energy in the vacuum takes the form: 

K(p) = S (0) + W- m)E« m) 2 S c (p) , (62) 

and being constant divergent quantities while S c (p) is a convergent 
function. Thus hard renormalization entails: 

5m = -E (0) , (63a) 
Z 2 -l = E (1) . (63b) 



Our scheme naturally separates the Hartree term (already discussed) and the 
Fock one: we get for the self-energy 

^v{p) — ^Hv + ^Fv(p) + convergent terms . (64) 



Similarly we rewrite the mass counterterm as 

5m = 5m H + 5m F + 0(h 2 ) , (65) 



3 We indicate the two- and more-points vertex functions in momentum space with 
a short-hand notation, which implies a Fourier transform. For example: 

,yi| 1 -(*, - i 2 y»<*'-»> — — — 



s (w/ d4( 



StpSi> (2vr) 4 J y ' 5ip(x 1 )5tl){x 2 ) 
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where the 0(h 2 ) terms are finite and 

5m H = -T, Hv ■ (66) 

Since Y>hv is independent of the momentum, it does not contribute to Z 2 — 1, at 
variance with S^„. The latter can be evaluated analytically but its complicated 
expression is not given since it is not used in the following. 

Next we come to the diagram c) of Fig. 3. Again, an explicit evaluation of this 
diagram is not required. The corresponding vertex function in the vacuum has 
the structure 

where is a logarithmically divergent constant and Vf\p, q) is a finite func- 
tion, vanishing at the renormalization point. Its choice is however a remarkable 
source of ambiguity, as we shall see for the a 3 vertex. In the present case we 
can safely choose to put all particles on-mass-shell. Furthermore, the explicit 
form of T^(p, q) is irrelevant in the "vacuum subtraction" scheme, since this 
amounts to neglect it at all. With the above definitions we get 

Z l - 1 = f (c) . (68) 



The previous considerations complete the renormalization program for the 
nucleonic part of the action. Now we consider the integration over the fermionic 
fields in the generating functional: 

Z[J] = ^J D[4>,ljj,a} e ifd 4 x[Cre n (x)+C c ,(x)+a(x)J(x)] _ (gg) 

This requires special care because it involves divergent quantities, which in 
addition are of different order in the BLE: we must keep Smn in the th order, 
being the only fermionic counterterm without mesonic loops, and leave the 
remaining counterterms inside the fermion loops that constitute the effective 
vertices of the bosonic action. 

To remove the explicit occurrence of Smn in the fermionic loops, the usual 
trick is to shift the scalar field as follows 

a^a = a-\ = a H -p sv (70) 

9 mi 

and replace in the total Lagrangian £ rcn + C ct a with a: 

C = Z 2 4>(ift — m)ip — SmpTpip + (Z 1 — l)5m H ijjip + Zigipaip 

+ ^ [d^a - ml& 2 ) - Umla 2 - ^Z,\^ - \z h \^ . (71) 
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We can now integrate over the fermionic fields to get the bosonic effective 
action 

- iTr log {-i [z 2 S Q l - 5m F + (Z 1 - l)5m H + Z ig <j] } . (72) 

As already shown in Section 3, the log in the previous action contains the 
closed fermionic loops, representing the effective vertices of our theory. They 
still contain, however, divergent counterterms, which need to be eliminated, 
as it is illustrated below. 

Along the same path of Section 3 we rewrite [using the symbol Aq of Eq. (25)] 

- iTr log {-i[Z 2 S 1 - 5m F + (Z 1 - l)8m H + Ziga}} 
= A* - iTr log {l + S [-5m F + (Z 2 - 1)S^ + (Z 1 - l)5m H + Z l9 a] } 

= A*+i ^Tr {^o [-Sm F + (Z 2 - 1)3^ + {Z 1 - l)5m H + Z l9 a] }" . 



n=l 71 



(73) 



Let us notice that —5m F + (Z 2 — l)^ 1 is the counterterm canceling the 
divergent part of the diagram b) in Fig. 3. Moreover the existence of the 
tadpole carries another contribution, [Z\ — l)5m,H, which compensates the 
divergence of the diagram of Fig. 4. 




Fig. 4. The tadpole with vertex correction. 



The last line of Eq. (73) requires some caution, since the appearance of the 
conterterms 5m F , Z 2 — 1 and (Z\ — l)5rriH seems to alter the order n of 
fermionic propagators So or, equivalently, of bosonic vertices. To clarify this 
point let us consider the Pauli-Villars regularization [4], which amounts to 
subtract from a diagram containing a a exchange (e.g. diagram b) of Fig. 3) 
a similar process with a heavy meson with mass A. The number of bosonic 
vertices remains obviously unaltered until one takes the limit A — > oo, which 
shrinks the insertion to one point only, corresponding to a counterterm. On 
this basis, as far as the number of bosonic vertices is concerned, one has to 
count twice the counterterm insertions of the type 5m F and Z 2 — 1 and, on 
the same foot, thrice the counterterms of the type Z\ — 1, because diagram c) 
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of Fig. 3 has 3 external points. This amounts to say that 5mp and Z2 — 1 are 
of order g 2 and Z\ — 1 is of order g 3 . 

The above considerations enable us to correctly count the number of bosonic 
vertices (namely of powers of g) in a fermionic loop like the one illustrated in 
the right of Fig. 5: indeed a loop with, say, n — 2 cr-vertices and a counterterm 
of the kind —5mp + (Z2 — l)^ 1 (left diagram in Fig. 5) is again of order g n . 
Within the "vacuum subtraction" scheme, the counterterms are the elemen- 




Fig. 5. Renormalization of the diagram b) insertions. 

tarily divergent diagrams evaluated in the vacuum; then the sum of the two 
diagrams of Fig. 5 reduces to the diagram on the right with the self-energy 
insertion regularized. The same occurs for the vertex corrections of diagram 
c) of Fig. 3 and of Fig. 4. 

The above procedure proves the renormalizability of the fermionic loops ap- 
pearing in the last line of Eq. (73) only for n > 5. Indeed for lower n's the 
fermionic loop itself is divergent and must be handled separately, by subtract- 
ing from these fermion loops their value at Jcf — 0. 

Let us then introduce the following notations: 

= 9 n (~i) tr / dxi . . . dx n ^"i^fe, • • • , x n )a(xi) . . . a(x n ) , (74) 



00 ^ j^n 

71=1 H 

A B V = E ^SS? , (76) 

n=5 n 

where the trace (the symbol tr) only refers to spin-isospin indeces. With the 
above definitions the bosonic action reads 

ATM far* ~ <° 2 ) ~ 

-h z ^ ~ h z ^} + E ^^0? +A* + A B M , (77) 
where the last two terms are now free from infinities. 
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o o - -6 



d) 



e) 



f) 



Fig. 6. The elementarily divergent diagrams contributing to the renormalization of 
d) the a meson mass and wave function and e) the 3a and f) 4a vertices. 

We then proceed to fix the renormalization parameters 5ml, Z 3 , Z 4 and Z 5 
in order to remove the remaining infinities. This goal is achieved in perturba- 
tion theory, starting from the meson self-energy diagram d) of Fig. 6, which 
renormalizes the a mass and coupling constant. The corresponding two-point 
vertex reads: 



5 2 T 



(78) 



where 



Ili d \ml) + (q 2 - ml 



dll^(q 2 

| v_12_ 

dq 2 



W-rnirn^iq 2 



(79) 



Il^(g 2 ) being a convergent function. Renormalization then requires: 

5ml = -U v (ml) 



Z.-l 



dli v (q 2 
dq 2 



q 2 =ml 



(80a) 
(80b) 



Let's consider in details the divergent contribution of diagram d) of Fig. 6. By 
applying Eqs. (80) we see that in Eq. (77) / d 4 x[^(Z 3 — 1) (d^ad^a — m 2 a a 2 ) — 
^5m 2 a a 2 \ and (1/2)2^ cancel each other up to finite terms which are vanishing 
on-shell as well as in the "vacuum subtraction" scheme. 



are given in the Appendix 



Explicit formulas for (ml) and — j^- 

together with their expressions both in the hard and in the soft renormalization 
scheme, the latter being adopted in Ref. [8], with renormalization point at 
q 2 = 0. As we shall see, the difference is relevant. 

In Fig. 6 the elementarily divergent diagrams contributing to the 3a [diagram 
e)] and 4a [diagram f)] vertices are also depicted. We define: 
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5 3 T 



<5 4 r 
da* 



(81) 
(82) 



where the quantities f ^ and are logarithmically divergent (in A) con- 
stants and r£ e ) and are suitable finite functions which disappear in the 
"vacuum subtraction" scheme. Renormalization enforces 



A 3 (^4 - 1) 
A 4 (Z 5 - 1) 



: f(<0 ^ 



and the divergent diagram e) [f )] of Fig. 6 exactly cancels the ZqJ [T^t 
tribution in the action (77). 



(83) 
(84) 

con- 



The above discussion does not exhaust the cancellation of infinities present 
in the action (77); other, non 1PI, divergent diagrams exist. This is linked to 
a subtlety in proving, via a Legendre transformation [5], that the functional 
T is the generator of the 1PI diagrams. A careful examination of this proof 
reveals that it does not apply to the tadpoles like T, Hv . For physical processes 
in the vacuum this does not create any pathology, since T, Hv is automatically 
renormalized to for hp = and no further divergences survive in (77). 
However, in the present context, the tadpole in the medium is neither vanishing 
nor trivial, and some care is required. Indeed we must also account for those 
diagrams that are not 1PI due to the presence of one or more T, Hv insertions. 
They contribute not only to F aa and T uaa but also to the term of the action 
which is linear in a and are shown in Figs. 7-9. We note that these diagrams 
are indeed originated by the quadratic, cubic and quartic terms in the shifted 
scalar field a = a — T<H V /g appearing in the effective action (77). When these 
diagrams are accounted for, taking care of their multiplicity, one sees that all 
the remaining divergences in the action cancel out. 



■O 



■o 



Fig. 7. Further divergent diagrams contributing to T c 



Let us consider in detail the terms of the action linear in a. The counterterms 
corresponding to the divergent diagrams of Fig. 9 cancel all these terms but 
one, namely — m 2 a a (^r 2 -)- The latter, however, is canceled by the term T^. 
Thus, the mean field solution in the vacuum is a(x) — 0, as it should. 
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Fig. 8. Further divergent diag rams contributing to r CT(T . 
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Fig. 9. Divergent diagrams contributing to the linear part in a of the bosonic effective 
action. 



Yet one consideration before leaving the problematic in the vacuum concerns 
the existence of many more divergent diagrams, which contain closed bosonic 
lines. However, the vacuum subtraction technique clearly renormalizes these 
contributions to zero. Finally, the renormalized bosonic effective action reads 



+ A B + A B M + finite terms . (85) 



The only difference with respect to Eq. (24) (besides the presence of the a 
self-couplings) is that ELi ^"[^off^]" is now replaced by constant, finite 
terms whose structure depends on the choice of the renormalization point. 



Before proceeding to study the self-consistency equation for the scalar field we 
observe that it is possible, as it was done in Ref. [8], to put these finite terms 
to by simply choosing the renormalization point at q 2 = 0. However, in so 
doing the a mass has to be re-evaluated order by order in the renormalization 
process. This would seriously plague a realistic theory (like, say, a model for 
the NN7r interaction) where the masses of the mesons as well as the scattering 
lengths are measured. In the context of QHD the mass of the a is not measured 
but plays the role of a phenomenological parameter which is fixed to reproduce 
the saturation properties of nuclear matter. Hence it is irrelevant whether the 
renormalization procedure shifts it, or not. There is however an important 
caveat: one could wonder if it has to be fixed before or after the renormalization 
process. We shall see that this point is of special relevance. 



We now derive the equation of motion for the mean field. Looking as usual to 
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a constant, uniform solution for this equation we get 
6A5r[a] 



5a 



= 

^3 2 ^43 . rp f d A q 



rrt^a - —a - —a 



2 6 y J {2-KYS,\q)+ga 

+ ^E(-l)VTr / JlS-S&Mo"- 1 = , 

n=\ J ( 27T ) 

where the last two terms in the lhs of the second equation are both diver- 
gent but their divergencies cancel each other. Within our cutoff regularization 
scheme: 

f d 4 q 1 

~ iTr / r9 x 4 Q -u \ 7 = Psijn- go) , (87) 

J (2tt) 4 S (q) + no- 
where the explicit expression of the scalar density p s is given in the Appendix. 

The remaining finite contribution in Eq. (86) can be evaluated by using the 
results of the Appendix and taking the limit A — > oo. By denoting with a the 
solution of Eq. (86) and recalling the effective mass definition (38) we get the 
self-consistency equation 

-mla - - ^a 3 + S sc (a) = , (88) 
2 6 

where 



£ sc (<x) = gp s {m -ga)+i £(-1) n g n Tr j ^-S^q)^' 1 = (89) 

9 \ ( -x3i k F + E* F 2 _ 5 _ 2 11 3 \ 

(m- gcr) log — (m - ga)k F E F -ma- -ma - —a' ^ 



and E* F = \Jk 2 F + (m - ga)^ 



Before concluding this discussion we wish to stress once more that the choice 
of the renormalization point introduces some ambiguity. Indeed it is an aspect 
of a more general question, namely the relevance that the quantum vacuum 
effects should, or should not, have on the model. 

In the previous derivation we have split the interaction brought in the bosonic 
action by the fermion loops into two parts: the density-dependent and vacuum 
contributions. We could reasonably ask whether the effect of the vacuum (the 
presence of the Dirac sea) should pertain to the particle or to the nuclear 
world. This matter is something more than a merely philosophical question, 
since it affects the self-consistency equation in a crucial way. If one indeed 
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adopts the vacuum subtraction scheme, one simply gets 

Ssc(o-) = gp S M(m - go) . (90) 

This eliminates the corrections (Relativistic Hartree Approximation) intro- 
duced in Eq. (89) and formally derived for the first time by Chin [8]. At the 
same time, the suggestion that the treatment of the vacuum dynamics in the 
approach of Ref. [8] is unsatisfactorily handled [32] seems to find a natural 
explanation. 

On the other hand, if one accepts the soft renormalization scheme of Ref. [8], at 
each perturbative order the a mass varies and its true value is reached when 
a fixed point, ruled by the renormalization group equation, is met. We can 
conclude that in the vacuum subtraction scheme, the a mass is fixed by hard 
renormalization; in this case the Relativistic Hartree contribution introduced 
by Chin must be neglected since it is an effect due to the dynamics in the 
vacuum. 



6 The Bosonic Loop Expansion 

In order to explicitly evaluate the partition function (43), which through 
Eqs. (18)-(22) supplies the static properties of nuclear matter, an approxima- 
tion scheme is required. As in Refs. [17,21], we use a Semiclassical expansion 
at the leading and next-to-leading orders: this amounts to evaluate first the 
partition function at the saddle point and then the quadratic quantum me- 
chanical fluctuations around the saddle point. The a self-interaction terms 
are disregarded in the formal derivations of the present Section. 

The standard Semiclassical approximation consists in an expansion in powers 
of H. By making explicit the dependence on this constant, one realizes that 
a factor TT l appears in front of the original action for fermions and bosons 
[see Eq. (11)]. However, after integration over the fermionic fields, Planck's 
constant also appears with a non-trivial dependence inside the new bosonic 
effective action of Eqs. (41) and (45). Thus, the use of the stationary phase ap- 
proximation does not lead to the standard Semiclassical expansion in powers 
of h. In this case, the ordinary way to proceed [5,6] is to introduce a dimen- 
sionless parameter gT 1 (a being small) in front of the bosonic effective action, 
to expand in powers of a and finally to set a = 1 at the end of the calculation 
at any given order in a. 

The bosonic effective action of Eq. (41) [equivalent to the Euclidean action 
(45)] has been constructed in such a way that the corresponding mean field 
equation only admits the solution a = 0: no term linear in a is present in 
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this action. Thus, we first separate in efr of Eq. (45) the terms quadratic in 
a from the higher order ones. By applying this procedure, from Eq. (43) we 
obtain the partition function in the form 4 

Z(/3)[J] = Z*{(3) e - n(S ^ 2 J D[a\ e 1 ai A ^ +A ^ + ^ I '**°Wm\ . (91) 

In the above, ° eS contains the part of the effective action which is quadratic 
in a and A^ eS the remaining effective interaction: 

00 / 1 \n 

AfM = ~ Tr S [She 9< ■ (92) 

n=3 n 

More explicitly, a first contribution to comes from the free bosonic field 

while a second one originates from the first term of the effective interaction in 
Eq. (45) and reads 

1 If 
--Tr[S HE a] 2 = --j d 4 xd 4 ya(x)tr[SHE(x~y)SHE(y-x)]a(y) (93) 

= -^Jd 4 xd 4 ya(x)UoE(x-y)a(y) , (94) 

which implicitly defines the equivalent of the Lindhard function in the rela- 
tivistic Euclidean case, n _E- The trace in the rhs of Eq. (93) (the symbol tr) 
only refers to spin-isospin indexes; in this case it just amounts to an implicit 
factor 4. Thus: 

AfA°\ = \° {Doe ~ U 0E ) a ee ^(A^V , (95) 

where we have introduced the Euclidean Random Phase Approximation (RPA) 
propagator 

Ai PA = (D J - iWf 1 , (96) 

with the diagrammatic representation given in Figure 10. Actually, only direct 
ring contributions enter into this propagator. 




Fig. 10. Dynamical content of the RPA-dressed boson propagator A RPA : dashed 
lines means free a propagators, solid circles represent fermionic loops. 



4 The introduction of the term with the external source J will be useful for later 
manipulations of the functional integral (see Eq. (101) and Section 6.3). The need 
for the factor yfa in front of the external source will become clear later. 
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6. 1 Mean Field 



The mean field contribution (er = 0) to the partition function (91) comes from 
the constant in front of the functional integral: 



,2 



It is well known that 



Z*((3)= J D$,$\ e J d A xd*yi,{x)S H l E (x-y)m (gg) 



p,A 



where E* = y p 2 + m* 2 , m* = m — ga and A is a spin-isospin index, is nothing 
but the partition function of a relativistic free Fermi gas up to the replacement 
of the bare nucleon mass with the effective one. The ground state energy of 
nuclear matter at the mean field level is then deduced from Eq. (18): 



E 



lim 



9 log Z(P)[J] 



mean field 



mean field fl— >oo Q(3 



+ EpA (99) 



777 

j2E;e(E F -E;) + n^a 

p,A 



(7 -2 



The well-known mean field result by Serot and Walecka [7] for the energy 
density: 



E 






mean field 



mean field 



l^^-^^ 2 (100) 



n J (2tt) 



k F E* F [(m*) 2 + 24] - (m*) 4 log ^ ' } + ^a 2 , 



An 2 [ m* j 2 

easily follows from Eq. (99) with the replacement J2p,\/& — > 4 / <i 3 p/(27r) 3 , 
the factor 4 being the spin-isospin degeneracy. 



6.2 One-loop correction 



In order to evaluate the functional integral in Eq. (91) we make the change of 
variable o — > \fa<J, thus obtaining 
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2 

*2 



Z(,8){J] = Z- W)e-"^'* (101) 

xvs/flwJ i ° (a|PAr, °- E -» 1 ^'' <? "' ,i[s " E! '° 1 "H 

2 

= Z*(/5)e-^ 2 {det[-Ai PA ]} 1 

where in the last equality a Gaussian integration has been performed 5 . 

The lowest order quantum mechanical correction beyond the mean field is thus 
obtained from the partition function: 

which, as we shall see, corresponds to the one-boson-loop approximation. By 
using Eq. (96), the one-loop contribution to the ground state energy of nuclear 
matter is obtained from: 



log Z((3)[J] 



(3 2f3 



if 00 1 

= — |Tr log[-Doi] - Tr £ -(D 0E U 0E r 

= -o I Tr hg[-D 0E ] - Tr / dX D ° eU ° E 



0E 



The first term in the last line is the vacuum energy of the a boson, which is 
washed out by renormalization. One also identifies the first term of the sum in 
the second line as the Fock contribution, in which the two nucleon propagators 
are Hartree-dressed (see Figure 1). The ground state energy at the one-loop 
level turns out to be 



5 Note that the Euclidean RPA propagator A PPA is negative-definite: 
Ai PA ( (? ) = -l/[ql + q 2 + ml + U 0E (q)] . 
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E 



lim 



log Z((3)[J] 



1— loop 



1-loop /3-^oo P 



(104) 



d 4 g } d ^ D 0E (q)u 0E (q) 



(2tt) 4 7 1 - XD 0E (q)U 0E (q) 



o 

where now, unlike Eq. (103), the fraction is truly algebraic. The above formula 
coincides (up to the isospin sum) with the Bruckner and Gell-Mann equation 
for the correlation energy of a degenerate electron gas in RPA approximation 
[33]. Moreover, it is also formally identical to the first-order approximation 
energy of the Modified Loop Expansion derived in Ref. [16] for a linear a 
model. It is evident that, up to the counting operator (namely the integration 
over A), Eq. (104) corresponds to a closed loop of the RPA-dressed boson 
propagator. 



6.3 Two-loop correction 



By expanding Eq. (101) up to first order in the parameter a, we find the 
following formal expression for the partition function at the next order in the 
Semiclassical expansion: 



Z(P)[J] 



2— loops 



l + ^T-/? (3) 



5J. 



(105) 



1(3) 



5 



J=0 



Here we have defined 



< # {3) (x 1 ,x 2 ,x 3 ) = S HE (x 1 - x 2 )S HE (x 2 - x 3 )S HE (x 3 - xi) , 



(106) 



^ 4 \x 1 ,x 2 ,x 3 ,x 4 ) = S H e(xi - x 2 )S H e{x 2 - x 3 )S H e{x 3 - x a )S H e{xa - Xi) , 
while «p (3) (^) 3 = Tr[S HE £j] 3 is a short notation for 

and an analogous definition holds for the term with four derivatives. 

The term proportional to ^Ja in Eq. (105) vanishes when setting the source J 
to and we finally obtain 



Z{(5)[J\ 



2— loops 



1 -a 



7 (3) 7 (3) 7 (3) 7 (4) 7 (4) 
— Z/ 1 — Zj 2 — Zj 3 -+- Zj 1 -T Zj 2 



(107) 
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where, analytically, 



(3) 



1 f 

z 3 3) = q9 6 J dx 1 dx 2 dx 3 ^ (3) {xi, x 2 , x 3 ) 

x J d Vl dy 2 dy 3 ^ 3) { yi , y 2 , y 3 )Al PA { Xl - Vl ) 

x|Al PA (x 2 -y 2 )Al PA (x3-y 3 ) 
+ Af A (x 2 -y 3 )Af A (x 3 -y 2 ) 
+3Af A (x 2 -x 3 )Af A ( y2 - y3 )\ , 



(108) 



and 



z[^ + = -g 4 J dx\ dx 2 dx 3 dx± (x\,x 2 , x 3 , X4) 
x{2Al PA ( a ; 1 -a; 2 )Al PA (x3-x 4 ) 



(109) 



+Af A (x 1 -x 3 )A^ A (x 2 -x,) 



RPA, 



In the last two equations, the different contributions to the functions zf^ and 
zf 1 appear, in the curly brackets, according to the order of the lhs terms. The 
Feynman diagrams corresponding to these functions are shown in Figures 11 
and 12. They clearly illustrate how the zf^ 1 and zf^ represent RPA-dressed 
two-boson-loop contributions. Note also that the Feynman diagrams of Fig- 
ure 11 (12) are the only distinct diagrams containing two boson-loops built 
up with three (two) boson propagators. 



Fig. 11. Feynman diagrams corresponding, respectively, to the contributions Z± , 
zf and Z {3) of Eqs. (108). 



At the two-loop level, the ground state energy of nuclear matter is immediately 
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Fig. 12. Feynman diagrams corresponding, respectively, to the contributions Z\ 
and of Eq. (109). 



(4) 



obtained from Eq. (107) by expanding log Z(f3)[J] 
and then setting a = 1: 



2— loop 



at first order in a 



E 



= — lim 

2— loops /3— +oo 

= lim — 



log Z((3)[J] 



2— loop 



(110) 



-z\ 



(3) 



,(3) 



zf + zr> + z. 



(4) 



,(4) 



Note that each one of the closed diagrams contributing to the zf^ and zf^ 
factorizes a (3VL. Indeed, due to translational invariance, the integrands in the 
various terms depend upon one less space-time variable than the number of 
space-time integrations. Therefore, the zero temperature limit in Eq. (110) is 
finite, as expected, and proportional to Q. 

Finally, from the two-boson-loop approximation for the partition function 



Z(P)[J}=Z((3)[J] 



mean field 



Z(P)[J] 



1— loop 



Z(P)[J] 



2— loops 



the ground state energy of nuclear matter is obtained as 



E = E 



mean field 



+E 



1— loop 



+E 



2— loops 



(111) 



(112) 



with contributions coming from Eqs. (99), (104) and (110). By applying Eqs. (20), 
(21) and (22), the binding energy per nucleon, pressure and compressibility of 
nuclear matter are thus derived. 



7 Numerical Results 



In order to illustrate the convergence properties of the proposed BLE, in this 
Section we present a few numerical results obtained in the one-boson-loop 
approximation. Despite the above formal development has been accomplished 
for a schematic dynamical model, the results discussed in the following corre- 
spond to a more realistic model, QHD-I, which includes the u meson beside 
the scalar-isoscalar one. In a forthcoming paper we shall discuss the formal 
modifications introduced by QHD-I into the schematic model considered here; 
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Table 1 

Coupling constants of the QHD-I model at the mean field and one-boson-loop 
approximations . 





9a 


9u 


Mean Field 


11.09 


13.81 


One-Boson-Loop 


13.00 


15.45 



we shall also systematically discuss the predictions for the static properties of 
nuclear matter of the a-cu model up to the two-boson-loop order. 

In our calculations we have adopted the following meson masses: m a = 550 
MeV and = 783 MeV. First we have investigated the properties of sym- 
metric nuclear matter obtained at the one-boson-loop level by adopting the 
values of the coupling constants g c and g u which allow us to reproduce the em- 
pirical saturation point of nuclear matter (B.E./A = —16.0 MeV and kp = 1.3 
fm -1 ) at the mean field level (see Table 1). For Fermi momenta larger than 
0.2 fm -1 the one-boson-loop binding energy turned out to be repulsive, with 
B.E./A = +18.6 MeV for kp = 1.3 fm -1 . This is mainly due to a strongly 
repulsive Fock term for the a which is not compensated by the attractive Fock 
contribution for the u. By adopting the same criteria for convergence estab- 
lished by Ref. [9] for the Loop Expansion and mentioned in the Introduction, 
we can conclude that the above large changes foreclose a "strong" convergence 
of the BLE. We then realise that even a non-perturbative expansion in the 
"large" couplings g„ and g u such as the BLE turns out not to be convergent 
in the "strong" sense. We remind the reader that the same finding was at- 
tained for the Loop Expansion [9], which, on the contrary, is perturbative in 
the couplings. 

It thus becomes important to investigate on a possible "weak" convergence 
of the BLE. To establish the validity of such a property, we have to check if 
it is possible to reproduce the empirical saturation of nuclear matter in the 
one-boson-loop approximation by a "reasonable" refit of the parameters. In 
Figure 13 we plot the binding energy per nucleon of symmetric nuclear mat- 
ter which we obtained in the one-boson-loop approximations by adopting the 
couplings of Table 1. The well known mean field predictions are displayed as 
well to facilitate the comparison with the result of the BLE at first order. Since 
the obtained refit modifies the couplings by a small amount (of about 15%), 
we can conclude that, contrary to the Loop Expansion of Ref. [9], the BLE is 
weakly convergent. Note that such an adjustment of the couplings practically 
maintains unaltered the underlying non-relativistic nucleon-nucleon interac- 
tion of the a-u model. 

We can therefore conclude that the scheme provided by the BLE proves to be 
a good candidate for performing systematic and reliable calculations, based 
on QHD models, beyond the tree-level approximation. 
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Fig. 13. Predictions of a a-to model for the binding energy per nucleon of symmetric 
nuclear matter obtained in the mean field and one-boson-loop approximations. The 
adopted coupling constants are given in Table 1. 



8 Conclusions and Perspectives 



In this paper we have proposed a functional integral approach to study the 
static properties of nuclear matter within a fully relativistic scheme, keeping in 
mind, even if not explicitly considered, the QHD model developed by Serot and 
Walecka. First, a bosonic effective action has been constructed for a system of 
interacting nucleons and a mesons. Special attention has been devoted to the 
formulation of a renormalization procedure for the introduced bosonic effective 
action. We have then evaluated the resulting partition function through a 
Semiclassical expansion around the solution given by the stationary phase 
approximation: this allowed us to express the result in terms of a boson loop 
expansion. Hartree-dressed nucleons and RPA-dressed bosons have revealed 
themselves to be the basic elements of this BLE. The suggested scheme has 
been analyzed at the mean field level and at the one- and two-boson-loop 
orders. We have presented a thorough exposition of the formal derivation of 
the binding energy per nucleon, the pressure and the compressibility of nuclear 
matter up to the second order in the BLE. 

Numerical results for the nuclear matter binding energy at the one-boson- 
loop level have been presented for a realistic, dynamical model including the 
oo meson in addition to the a. The convergence properties of the BLE have 
been discussed in this context. A systematic calculation of the static properties 
of nuclear matter up to the two-boson-loop order and in the a-u model will 
be the object of a forthcoming paper, where we shall also deal with the formal 
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derivation of the higher orders of the BLE. The inclusion of the u meson is not 
a prohibitive task from the point of view of the diagrammatic content, since 
the basic topology of the boson loops (see Figures 11 and 12) is common to 
any meson field which can contribute to the dynamics. However, since the a; 
component is coupled, in infinite nuclear matter, to the a meson (due to the 
breaking of the Lorentz invariance induced by the infinite medium) [34], the 
RPA equations, which deeply intervene in the corrections beyond the mean 
field, are much more involved that those presented here. 

Further, we note that the implementation of the isovector mesons p and 
5 [/o(980)] in our framework does not introduce additional difficulties from 
a technical point of view. Up to now, however, they have been handled only 
at the mean field level [27]. 

It is worth noticing that at the mean field level no pseudoscalar and pseu- 
dovector mesons can be coupled to the nucleon density or current, since the 
average values of pseudoscalar and pseudovector quantities vanish. However, 
the situation is completely different when boson loops are involved. It is indeed 
known that the two-pion exchange diagrams, with the possible simultaneous 
excitation of one or two intermediate nucleons to a A, provide a large amount 
of the nuclear binding. We shall explore this mechanism in the future, but we 
reasonably expect that a drastic change in the underlying dynamics of QHD 
will occur and that the gap between QHD and Bonn potential models will be 
further, and significantly, reduced. 



Appendix 



In the regularization scheme expressed by Eq. (59) the explicit expression of 
the scalar density, p s (m) = p sv {m) + p S M(m), is given by: 



. , mA 4 f m 2 A 2 1 

Psv{m) = . t^t \ 1 + — log — \ , 

l-n 2 [m 2 — A 2 ) [ m 2 — A 2 m z J 
PsM{m) = < k F E F - m log 



7T 2 (m 2 — A 2 ) 2 [ /;» 

with the limiting behavior 

m ( . o „ o , A 



m / A \ 

Psvijn) — > -— A 2 - 2m 2 log — + m 2 . 

A^oo 2n 2 V m J 



The quantity 11^ (q 2 ) introduced in Eqs. (78) and (79) reads 



37 



2ir 2 q(m — A) 2 (m + A) 2 
-(-m 6 + llg 2 m 4 + 3A 2 m 4 - 7g 4 m 2 - 3A 4 m 2 + 6g 2 A 2 

+m 2 + g 6 + A 6 _ g 2 A 4 _ g 4 A 2^ 

A 4 



x 



X 



An 2 q 2 {m 2 — A 2 ) 2 y/m — q — Ay/m + q — Ay/m — g + Ay/m + g + A 

2 2 

m z — g z 



arctanh 



,2 _ „2 _ A 2 



-arctanh 



V m — g — Ay/m + q — A\/m — q + A^m + g + A 

m 2 + g 2 - A 2 \ 
^ a/?ti — g — Ay/m + g — Ay/m — g + Ay/m + g + A y 



(m 4 - 6g 2 m 2 - 2A 2 m 2 + g 4 + A 4 ) log (^) A 4 



47r 2 g 2 (m 2 - A 2 ) 2 



47r 2 (m 2 - A 2 ) 



In the limit A — > oo and within a hard normalization scheme we find 



nf(m 2 ) 



- i -JA 2 -(6m 2 -m 2 )log- 
Z7r 2 I m 

4m 2 + m 2 (4m 2 — m 2 )i 
H - + — arctan 



^4m 2 



mi 



and 



ffl( rf )(g 2 ) 
dq 2 



1 2m 2 



2tt 2 I 4 
J 4m? - m 2 (2m 2 + m 2 ) 



arctan 



mi 



^4 m 2 _ m 2 



log 



On the contrary, by choosing the renormalization point at g 2 = we find 

TTMWnN A4 f 3m 2 + A 2 o 2 3A 2 + m 2 , Al 
v 2ix 2 [(A 2 -m 2 ) 2 (A 2 -m 2 ) 3 m) 

— > ( A 2 - 6m 2 log — + 5m 2 1 , 

A->oo 27T 2 I m J 
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and 



dnj d) (q 2 ) 
dq 2 



-A 4 (13A 4 - 68m 2 A 2 - 41m 4 ) 



g2=0 | 24tt 2 (A 2 - m 2 ) 4 

A 4 (A 6 - 3m 2 A 4 - 12m 4 A 2 



+ 



27r 2 (A 2 -m 2 ) 5 



2m 6 , A 
log — 

m 



13 1 i A 

A^ _ 24^r 2 + 2^° g m 



For the aaa vertex function of Eq. (81) we have 

2mA 4 (m 2 + 2A 2 ) . mA 4 (m 4 + 8A 2 m 2 + 3A 4 ) 



f(e) 



7r 2 (m - A) 3 (m + A) 3 
4m 3m , A 



+ 



7r 2 (m- A) 4 (m + A) 4 

; 3m A 
^ + — log — 

and for the a aaa self-coupling of Eq. (82): 
(7m 4 + 34A 2 m 2 + 7A 4 ) A 4 



log 



A 



m 



f(/) 



3vr 2 (A -m) 4 (m + A) z 

7 1 , A 
+ — !og — • 

A^oo Stt z 7f z m 
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